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We study the effect of the Singwi-Tosi-Land-Sj Slander local field correction on spin relax- 
ation/dephasing in n-type GaAs quantum wells at low temperature by constructing and numer- 
ically solving the kinetic spin Bloch equations. We calculate the local field factor G(q) in quantum 
wells by numerically solving three equations which link the local field factor, the structure factor, 
and the dielectric function, self-consistently. Such a correction reduces both the electron-electron 
Coulomb scattering and the Coulomb Hartree-Fock term. We compare the spin relaxation times 
with and without this correction under different conditions such as temperature, electron density, 
well width and spin polarization. We find that this correction leads to a decrease/increase of the 
spin relaxation time in the strong/weak scattering limit. At high spin polarization, it reduces the 
Hartree-Fock term and consequently tends to decrease the spin relaxation time. The modification 
of the spin relaxation time by the local field correction is more or less moderate either due to the 
coexistence of scattering other than the Coulomb scattering at low spin polarization and/or due 
to the competing effects from the Coulomb scattering and the Coulomb Hartree-Fock term at high 
polarization. 



PACS numbers: 72.25.Rb, 71.10.-w, 67.30.hj, 77.22.-d 
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The study of spin relaxation/dephasing (R/D) in semi- 
conductors is of great interest to the spintronic soci- 
ety recentlyfi^ It is well known that the D'yakonov- 
Perel' (DP) mechanism 4 is the dominant spin R/D mech- 
anism in both bulk and two-dimensional n-type zinc- 
blende semiconductors.— This mechanism is due to the 
momentum-dependent spin state splitting, which origi- 
nates from the Dresselhaus spin-orbit coupling^ in crys- 
tals without inverse symmetry and the Rashba spin-orbit 
coupling 7 - in quantum wells (QWs) with asymmetric po- 
tential. For (001) GaAs QW with small well width, the 
Dresselhaus term can be written as^ 
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in which (Arj) represents the average of the opera- 
tor — (d/dz) 2 over the electronic state of the lowest 
subbandi^ 

It was first pointed out by Wu et al. from a fully mi- 
croscopic approach that any type of scattering, including 
the spin conserving electron-electron Coulomb scattering, 
can lead to irreversible spin R/D in the presence of in- 
homogeneous broadening ] 10 : 11 ! 12 : 13 ! 14 Similar claim was 
also made later by Glazov and Ivchenkoi 1 ^ This has been 
verified experimentally by Ley land et al.. 16 Moreover, 
Weng and Wu also predicted an interesting effect that 
the spin R/D can be suppressed by increasing the initial 
spin polarization due to the effective magnetic field which 
originates from the Hartree-Fock (HF) contribution of 
the Coulomb interaction.— This prediction has also been 
confirmed experimentally very recently. 17 However, in all 
the previous works ) 9 ' 10 ! 11 ' 12 : 13 ' 14 ! 17 the electron-electron 



Coulomb interaction is treated under the random phase 
approximation (RPA). Usually the RPA is considered 



to be good only at small 
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with r, = 



representing the dimcnsionlcss coupling parameter for 
the two-dimensional electron gas with electron density 
n and effective Bohr radius a B . Such kind of approxima- 
tion does not involve the correction which describes the 
electron charge depletion nearby the electron, known as 
the exchange-correlation hole. This correction beyond 
the RPA is so-called the local field correction (LFC). 
It must be taken into account for properties when the 
electron-electron interactions plays an important role, 
eg., the spin Coulomb drag effect fi 9 - This LFC can be 
introduced through a momentum-dependent correction 
factor G(q) in the dielectric function^^^ 2 ^ 4 * 2 ^* 2 ! 
e(q,w) = l-« g pW(q,w)/[l + w g G( ? )p( 1 )(q,cj)], where 
P^ 1 ' (q, oj) is the free-electron dynamical polarizability 
and v q = 47re 2 /q 2 is the bare Coulomb potential. This 
factor was originally introduced by Hubbard 20 - in very 



simple form G H (q) = | q 2 q +k -i with fcj? being the Fermi 
wave vector. This form is unappreciated nowadays due 
to the lack of accuracy. Later, several improved forms 
of G(q) based on the diagrammatic techniques were 
proposed i 21 : 22 Singwi et al£^ gave a different and more 
precise choice of local field factor G STLS (g) as a functional 
of static structure factor 5(g) by using the equation of 
motion method: 
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They suggested that one can solve a triad of equations 
which link the three quantities: G STLS (g), the structure 
factor 5(g) and e(q, to) self-consistently. Later on, there 
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are many other efforts on calculating this factor in three- 
dimensional (3D) systemsi2&>2£>2M£i2£ii2iiiS In ideal two- 
dimensional (2D) electron gas/liquid, the local field fac- 
tors G 2D (q) were also calculated^^>^i2£i21i2&2£ similar 
to the 3D case mentioned above. 

In this paper, we calculate G w (q) in n-type GaAs QWs 
by solving three equations which link the three func- 
tions G w (q), S w {q), ande w (q, w) self-consistently. These 
equations which are a little different from those in the 
bulk system 2 ^ and ideal 2D system^ can be written as: 



e w (q,c) = l 
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where v q = V Q \I{i<lz ) | 2 and vq — 4ne 2 / Q 2 with Q = 
(q, q z ) are the screened quasi-2D and bare bulk Coulomb 
potential respectively. \I(iq z )\ 2 denotes the form factor 
of square QWs with finite well depth, whose expression 
can be found in Ref. 

We plot the wave- vector dependence of G w (q) in 
GaAs/Alo.4Gao.6As QWs with well width a = 10 nm 
under different conditions in Fig. [TJ In Fig.QJa), G w (q) 
are plotted at different temperatures (T) when the elec- 
tron density n = 2.5uq — 2.5 x 10 11 cm -2 (the corre- 
sponding r s = 1.09). One can see that G w (g) increases 
very quickly with temperature. The electron-density de- 
pendence of G w (g) is plotted in Fig. QJb) at T = 30 
K, where the largest electron density is chosen to be 
5no (r s — 0.77), still beyond the RPA limit. It is seen 
that G w (q) decreases with electron density n. We fur- 
ther present G w (q) at different spin polarizations P when 
n = 2.5no and T = 30 K. One finds that G w (g) increases 
slightly with P. When we use the infinite-well-depth as- 
sumption and let the well width to be very small, our 
G w (g) tends to G 2D (q)M- 

Facilitated with the local field factor, we turn to in- 
vestigate the effect of LFC on spin dephasing in n-type 
(001) GaAs/Alo. 4 Gao. 6 As QWs. We construct the ki- 
netic spin Bloch equations by using the nonequilibrium 
Green function method: 14 
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where pu^a' denote the single particle density matrix el- 
ements with the diagonal and off-diagonal elements to 
be the electron distribution functions /k CT and spin co- 
herence pk,<T-<r- The coherent terms pk,<r<r'\coh describe 
the precession of the electron spin due to the effec- 
tive magnetic field fi(k) [Eqs. <$$5§] as well as the HF 
Coulomb interaction^ In the scattering term f>k,acr' \scatu 
the electron-LO-phonon, the electron-AC-phonon, the 
electron-nonmagnetic impurity and the electron-electron 
Coulomb scattering are included explicitly. Their expres- 
sions can be found in Refs. 0[T2|. After numerically solv- 
ing the kinetic spin Bloch equations, one can obtain the 




FIG. 1: (Color online) The local field factor G w (q) versus the 
wave vector in GaAs/Alo.4Gao.6 As QWs with well width a = 
10 nm. (a) At different temperatures when electron density 
n — 2.5no; (b) For different electron densities at T = 30 K; (c) 
For different spin polarizations at T — 30 K and n — 2.5no- 



spin dephasing and relaxation times from the temporal 
evolutions of the spin coherence p^ a -a and the electron 
distribution function /k )( r>^ 

In Fig. [2j we present the temperature dependence of 
the spin relaxation time (SRT) in QW with (t L fc) and 
without (t) the LFC for two different electron densi- 
ties n = 2.5no and 5hq. The spin polarization is small 
(P = 2%) in the calculation. The well width a — 10 nm. 
We choose the impurity density to be 0.0084n , corre- 
sponding to that we used in Ref. Il7l . In our calculation, 
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FIG. 2: (Color online) SRT versus temperature with (solid 
curves) and without (dashed curves) the LFC for two different 
electron densities n = 2.5no and 5no. a = 10 nm and P = 2 
%. 



we focus on low temperature case. It is seen from the fig- 
ure that t lfc is a little larger than r at very low tempera- 
ture and becomes smaller than r at higher temperatures 
for both cases. These results are due to the reduction of 
the Coulomb scattering. It has been pointed out that the 
SRT decreases with the total scattering strength (1/r* 
with r* representing the momentum relaxation time^^) 
in the weak scattering limit (|0|t* > 1), but increases 
with the total scattering strength in the strong scattering 
limit {\U\t* < l). 13 »tH Therefore, at very low tempera- 
ture, the system is in the weak scattering limit and the 
decrease of the Coulomb scattering leads to t lfc > r. At 
higher temperature, however, the system is in the strong 
scattering limit. Consequently, with the decrease of the 
Coulomb scattering, t lfc < t. It should be noted that 
the amplitude of the modification of the SRT is no more 
than one order of magnitude because of the existence of 
other scattering terms such as the electron-phonon and 
the electron-impurity scattering. 

From Fig. [21 one also finds that the reduction of the 
SRT by the LFC for the case with n = 2.5uq is much 
larger than that with n = 5uq. Therefore, it is necessary 
to further investigate the density dependence of the SRT 
with/without the LFC. In Fig. the SRTs are plotted 
against the electron density n for different well widths 
a = 10 nm and 20 nm when T = 30 K and P = 2 %. 
In order to make the figure clear, we re-scale the SRTs 
for the case with a = 20 nm by multiplying 0.2. As the 
system is in the strong scattering limit, r LFC < r. More- 
over, the larger the electron density is, the smaller the 
difference between r LFC and r becomes. This is because 
the reduction of the Coulomb scattering due to the LFC 
becomes weaker for larger electron density. This can be 
seen from Fig.QJb) that G(q) decreases with n, especially 
at large q near the Fermi momentum fc^- 

Finally, we investigate the spin polarization depen- 
dence of t and t lfc for different temperatures and elec- 
tron densities as shown in Fig. In the calculation, the 



FIG. 3: (Color online) SRT versus electron density with (solid 
curves) and without (dashed curves) the LFC for two different 
well widths a = 10 nm and 20 nm. T = 30 K and P = 2 
%. The SRTs for the case with a = 20 nm are re-scaled by 
multiplying 0.2. 
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FIG. 4: (Color online) SRT versus spin polarization with 
(solid curves) and without (dashed curves) the LFC for two 
temperatures T = 10 K and 30 K and two electron densities 
n = 2.5no and 5no. a — 20 nm. 

well width a = 20 nm. It is interesting to see that there 
is a crossover at certain spin polarization. r LFC < T for 
small spin polarization and r LFC > r for large spin polar- 
ization. These features origin from the dual effects of the 
Coulomb interaction at high spin polarization. The first 
effect has been said above. The second one is from the 
effective magnetic field in the HF term which suppresses 
the spin R/D at high spin polarization. — This effect has 
been observed in experiments very recentlyj 17 i 41 As the 
LFC always reduces the Coulomb interaction, the LFC 
can lead to a decrease of the SRT by reducing the HF 
term at high spin polarization. It is seen from the figure 
that for small spin polarization, as the contribution from 
the HF term is negligible and the system is in the strong 
scattering limit, inclusion of the LFC always results in a 
decrease of the SRT. However, when the polarization be- 
comes large enough, the contribution from the HF term 
becomes important and the system turns into the weak 
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scattering limit as the HF term enhances the amplitude 
of effective magnetic field |fi| along the z-axisi 11 ' 17 In 
this regime, inclusion of the LFC results in a compet- 
ing effects from the HF term and the scattering term. 
Our results indicate that when the polarization is large 
enough, the effect from the scattering is stronger and 
consequently the SRT is increased with the inclusion of 
the LFC. 

In conclusion, we study the STLS LFC on the spin 
relaxation/dephasing in n-type GaAs (001) QWs at low 
temperature by constructing and numerically solving the 
kinetic spin Bloch equations. The LFC takes into account 
of the screening from the additional exchange-correlation 
hole ignored in the RPA. The LFC reduces both the 
Coulomb scattering and the Coulomb HF interaction. 
It is known that the decrease of the Coulomb scatter- 
ing can lead to a decrease (an increase) of the SRT in 
the strong (weak) scattering limit whereas the decrease 
of the Coulomb HF interaction results in a decrease of 
the SRT at high spin polarization. 11 After comparing the 



SRT with and without the LFC under different condi- 
tions, we find that the LFC can lead to a decrease of the 
SRT in the strong scattering limit for small spin polar- 
ization. It also leads to a decrease of the SRT while the 
spin polarization getting larger, and finally leads to an 
increase when the spin polarization is large enough. 

It is noted that in our recent theoretical comparison 
with the experiment, the LFC was ignored^ 7 - From this 
investigation, we conclude that this approximation is ac- 
ceptable. This is because when the electron density is 
about 2.1no, the difference between r LFC and r is very 
small according to Fig. [4] 
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